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ABSTRACT 

The extragalactic background light (EBL), i.e., the diffuse meta-galactic photon field in the 
ultraviolet to infrared, is dominated by the emission from stars in galaxies. It is, therefore, 
intimately connected with the integrated star formation rate density (SFRD). In this paper, 
the SFRD is constrained using recent limits on the EBL density derived from observations 
of distant sources of high and very-high energy gamma-rays. The stellar EBL contribution is 
modeled utilizing simple stellar population spectra including dust attenuation and emission. 
For modeling the SFRD up to z = 4 a broken power law function in z -H 1 is assumed. A wide 
range of values for the different model parameters (SFRD(z), metallicity, dust absorption) is 
investigated and their impact on the resulting EBL is studied. The calculated EBL densities 
are compared with the specific EBL density limits and constraints on the SFRD are derived. 
For the fiducial model, adopting a Chabrier initial mass function (IMF) and a second power 
law index for the SFRD of (5 - 0.3, the SFRD is constrained to < 0.1 Moyr ' Mpc"^ and 
< 0.2 Mq yr"' Mpc"^ for a redshift of z ~ 1 and z ~ 2, respectively. The limits for a redshift 
of z ~ 1 are in tension with SFRD measurements derived from instantaneous star formation 
tracers. While the tension for the conservative fiducial model in this study is not yet overly 
strong, the tension increases when applying plausible changes to the model parameters, e.g., 
using a Salpeter instead of a Chabrier IMF or a adopting a sub-solar metallicity. 
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1 INTRODUCTION 

The star formation rate density (SFRD) describes the evolution of 
stellar formation over the history of the universe and is closely con- 
nected to structure formation and reioniziation (for a recent review 
see [Robertson et aLl|2010^ . The SFRD can be obtained by com- 
bining the star formation rate for individual galaxies as derived, 
e.g., from the total infrared (IR), radio luminosity, or the Lyman-or 
emission (Kennicutt 1998|, with the space density of the sources 
(see |Hopkins,2004,|Hopkins & Beacom|2006| for recent data com- 
pilations and references). 

The extragalactic background light (EBL; see |Hauser & Dwek] 
|200H for an excellent review), which is the diffuse meta-galactic 
photon field in the ultraviolet (UV) to the infrared (IR) wavelength 
regime, is dominated by stellar emission in the optical (O) to near- 
infrared (NIR) and by stellar emission reprocessed by dust in the 
mid to far-infrared (MIR/FIR). It, therefore, provides a probe of 
the integrated SFRD. Direct measurements of the EBL are chal- 
lenging due to strong foreground emission in our planetary system 
(zodiacal light) and galaxy ([Hauser et al.|1998]). Lower limits on its 



density are derived from integrated source counts (e.g. [Madau &| 
[Pozzet ti 2000; F azio et al.|2004[|Dole et al.|2006> . 

The most constraining upper limits on the EBL density are ob- 
tained from spectroscopic observations of distant sources of very- 
high energy gamma-rays (VHE; £ > 100 GeXTlc VHE gamma-rays 
produce electron-positron pairs with the low energy photons of the 
EBL (7vHE 7EBL — > ^ e*), resulting in an energy dependent attenu- 
ation signature in the measured VHE gamma-ray spectra ( |Nikishov| 
[1962) |Jelley|[T966l [Gould & Schreder|[T967) . Using assumptions 
about the intrinsic spectrum emitted at the source the EBL density 
can be constrained (Stecker et al.|199"2l ). The discovery of distant 
(z ~ 0.2 - 0.5) sources of VHE gamma-rays with current gener- 
ation ground-based VHE instruments like H.E.S.S., MAGIC, and 
VERITAS lead to strong constraints on the EBL density, in partic- 
ular in the O to NIR (AharoniaiTet al7 20061 [Albert et al.[[2008l >. 
These instruments also increased the number of known extragalac- 
tic VHE sources from about five in the year 2000 to more then 50 
in the year 201 Ip] By combining the observations of several VHE 
sources, the systematic uncertainties in deriving upper limits on the 
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EBL density can be significantly decreased and constraints over a 
wide wavelength range of the EBL can be derived (Dwek & Kren-| 
nrich 2005 Mazin & Raue|2007} . In particular, in a recent study by 
Meyer et al. p012 ME 12 in the following) data from the Fermi- 
LAT instrument, covering the MeV to GeV energy regime, were 
combined with a large sample of extragalactic VHE sources to de- 
rive the strongest and most robust upper limits on the EBL density 
today. 

Limits and measurements of the EBL can be used to investi- 



gate the SFRD in the early universe ( [Santos et al. 2002, Fernandez 
& Komat su'2006, Raue et al. 2009! |Maurer et al.|2012||Gilmore 



2012| >. In principle such studies provide a way of investigating the 
aggregated star formation of galaxies too faint to detect individu- 
ally with current technology. Unfortunately, the current limits are 
not strongly constraining, due to the weak expected contribution 
to the EBL from such sources and the lack of sufficiently precise 
measurements of the EBL density. 

The dominant contribution to the EBL density comes from 
star formation at redshifts of z ~ 1, where a peak in the SFRD is 
expected (e.g. |Hopkins & Beacom|20^6l l. The EBL is, therefore, an 
excellent probe of the bulk star formation (e.g. |Dwek et al.|1998 
Madau & Pozzetti|2000[|Chary & Elbaz|2001[ >. For example, |FMdall 
et al.| P007J studied the SFRD by combining the integrated EBL 



density with measurements of the observed stellar mass density. 
They concluded that the measurements for these two quantities are 
in tension, which could be weakened by adopting a different univer- 
sal stellar initial mass function (IMF). [Horiuchi et al.| ( |2009t used 
the integrated EBL density calculated for specific SFRDs together 
with measurements of the EBL to constrain their modeling of the 
diffuse neutrino background from supemovae. 

In this paper, the SFRD will be investigated utilizing limits on 
the EBL density derived in ME12. In particular, conservative upper 
limits on the SFRD will be derived. To this end, the EBL density 
is calculated using emission spectra from simple stellar population 
(SSP) modeling (Section [2|. The employed methodology is simi- 
lar to the one presented in Raue et ^ ( [2009^ , but it is extended by 
making use of the full wavelength dependent EBL constraints and 
includes modeling of the MIR/FIR emission from dust. The EBL 
model is constructed to be in minimum tension with the upper lim- 
its to derive conservative constraints on the SFRD. The effect of the 
different parameters entering the model are investigated and upper 
limits on the SFRD are derived (SectionHb. The implications of the 
the results are discussed in Section]?] 

Throughout this work a flat ACDM cosmology with £1^ = 0.7, 
Q.U = 0.3, andZ/o = 70km s"' Mpc"' ('737') is assumed. 
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Figure 1. Resulting EBL SED (z = 0) for different choices of parameters of 
tlie EBL model. Upper Panel: IMF and SSP models. Lower Panel: Metal- 
licities. To guide the eye, various EBL measurements and upper and lower 
limits are shown as grey markers in the background (from |Hauser & Dwek| 
laOOTllMazTn & Raue|2007||Raue & Mazin|2011> . 



yields the energy density, or, after multiplication with vcjAn, the 
co-moving spectral energy distribution (SED) of the EBL 



Py{z) = ViM = V 



An 



i: 



GAz') 



df 
dt 



dz' 



(2) 



with v' = y(l +z')/(l +z). 

In the following, the various model parameters are discussed. 



2 EBL MODEL CALCULATIONS 



The EBL density is calculated following |Dwek et al.| ( |1998] l and 
[Kneiske et"ar| (J2002) (KN02 in the following) from the star for- 
mation rate density p, (SFRD) and the specific luminosity Lv(t) of 
a simple stellar population (SSP) model of age t as, e.g., derived 
from stellar population syntheses models. The co-moving emissiv- 
ity (luminosity density) at redshift z is obtained from the convolu- 
tion 



fiv(z) ■■ 



r 



LMz) - t(z'))p,(z') 



dz' 



(1) 



where the SFRD p,(z) is assumed to start at some finite epoch z,„ = 
z(t„i) and t(z) is the cosmic time corresponding to a redshift z- For 
given evolution of the emissivity a second integration over redshift 



2.1 IMF, metallicity, and SSP model 



The SSP spectra used for the EBL modeling are taken from |Bruzual| 
|& Cha riot ( 2003 1 (BC03 in the following). 

IMF The initial mass function (IMF) is one of the fundamental pa- 
rameters defining the emission properties of the SSP spectra (see, 
e.g., [Kroupa et al.|20lT| for a recent review on the stellar IMF). 
Experimentally determining the IMF is a challenging tasks. A uni- 
versal IMF has been found from studies of stellar populations in 
the Milky Way (Kroupa IMF, |Kroupa|200f| , but recent studies of 
other galaxies indicate that low-mass stars might be more numer- 
ous in late-type galaxies (van D okkum & Conroy|20To| and that in 
early-type galaxies the IMF varies systematically with the mass-to- 
light ratio ^Treu et al.|2010[ [Cappellari et al.|2012|). In particular. 
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|Cappellari et al.| ^2012^ found that the IMF for a large sample of 
early tjfpe galaxies vary from a Kroupa to a Salpeter IMF lSalpeter| 
|1955[ l, providing more low mass stars (m < 0.5 Mq). 

BC03 provide SSP model spectra for a Salpeter (Salp) and 
a Chabrier (Chab) IMF (an IMF very similar to the Kroupa IMF; 
for a comparison see, e.g., Figure 28 of [Kroupa et al. 2011^ in the 
mass range from 0.1 to 100 Mq. The resulting EBL densities from 
the model calculations for a Salp and a Chab IMF are compared in 
Figure [T] Top Panel. The EBL density for a Chab IMF is ~ 30 % 
higher then the one resulting from a Salp IMF (all other model 
parameters are kept the same). This is an effect of the break of the 
Chab IMF at ~ 0.5 Mq: for a mass normalized IMF a Chab IMF 
produces more stars in the region above 1 Mq compared to the Salp 
IMF, which does not show this break (see Figure 28 of [Kroupa] 
|et al.|[20ri I. Since the bulk of the emission is produced by stars 
with > 1 Mq adopting a Chab IMF results in more luminosity per 
mass. 

A Chab IMF will be adopted for the fiducial model, but results 
for a Salp IMF will also be presented. While a Chab IMF will result 
in stronger absolute constraints on the SFRD it turns out that, due 
to the different behavior of instantaneous and late-time tracers for 
the SFRD for a specific IMF (see Appendix [A|, adopting a Chab 
IMF results in weaker constraints relative to the SFRD derived from 
instantaneous tracers. 

Metallicity The history of the metallicity of stars and star forming 



regions can be investigated using stars in our galaxy (e.g. Edvards 
son et al^|1993^ or large samples of distant galaxies (e.g. Panter 



et al. 20081. These studies show that the solar metallicity value 



of Zq = 0.02 is close to the average metallicity over time up to 
redshifts of at least z = 2. The metallicity evolution is a function 
of galaxy mass: while massive galaxies (M > 10" Mq) are en- 
riched early and show an almost constant metallicity since redshift 
of z ~ 2, lower mass galaxies show a decrease in metallicity l |Pan-| 
|ter et al.|2008^ . This decrease in metallicity can also be seen in the 
stars of our galaxy, where old stellar populations show lower metal- 
licities (e.g. |Edvardsson et al.|I993] l while slightly higher metallic- 
ities then the solar metallicity (up to [Fe/H]~ 0.2) are observed in 
present day star forming clouds ( [Binney & Merrifield|1998[ l. Over- 
all, the average metallicity at redshifts greater than 0.2 is well de- 
scribed by a solar metallicity or less ( [Panter et al.|2008} . 

The impact of different metallicities on the resulting EBL den- 
sity is displayed in Figure [T| bottom panel. A sub-solar metallic- 
ity value leads to an increased emission in the UV to NIR up to 
wavelengths of ~ 1 - 2/jm due to the lower mass loss in massive 
stars and hence increased time-integrated emissivity. A metallicity 
higher then the solar value leads to an increased EBL density in the 
NIR but less emission in the UV to O. A solar metallicity is adopted 
as the fiducial value for the EBL modeling and the effect of lower 
and higher metallicities on the resulting constraints on the SFRD 
will be investigated. 

SSP model For comparison. Figure [Tl Top Panel also shows the 
resulting EBL SED adopting SSP spectra calculated with the Star- 
burst99 code (Version 6.0.3;'Leithe rer et al.| 1999 [[Vazquez & Lei-[ 
[therer|2005] [Leitherer et"a l. 2010), for a Salp and a Kroupa IMF. 
Excellent agreement is found with the calculations using the BC03 
SSP spectra, with only small difference in the UV. These could arise 
from the different treatment of, e.g., the thermally pulsing asymp- 
totic giant-branch (TP-AGB) phase (see, e.g., [Bruzual & Charlot[ 
[2003i for a more detailed comparison of the models). An increase 
of the overall luminosity in the TP-AGB phase as, e.g., predicted by 
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Figure 2. Resuhing EBL SED (z = 0) for different choices of the infrared 
conversion factor /,- (upper panel) and the color of the dust attenuation 
E(B- V) (lower panel) in comparison the EBL upper limits from ME12 (red 
thick line). Shown are results for a Chab IMF, £(B-V) = 0.15,/;,- = 3x10', 
and a SFRD with /? = and zo = 1 if not denoted otherwise in the legend. 
To guide the eye, various EBL measurements and upper and lower limits 
are shown as grey markers in the background (from [Hauser & Dwek|200T| 
jMazin & Raue|2007|[Raue & Mazin|201 1\ . 



[Marastoii| ( [2005[ l, would increase the resulting EBL and strengthen 
the limits. 



2.2 Dust absorption and emission 

Part of the stellar emission is reprocessed by dust. Dust attenuation 
is calculated following KN02 and applied directly to the SSP spec- 
tra. Full absorption of all ionizing photons is assumed with 50?^of 
the absorbed photons being reemitted as Lyman-o- emission, which 
is consequently absorbed by dust. In addition, an average extinction 
curve is applied: 



Aa = 0.68 ■E{B-V)-R- (^"' - 0.35) 



(3) 



with R = 3.2, A being the wavelength in /jm, and the absorption 
coefficient being given by g(A) = lO""'''*'' . KN02 assumed two dif- 
ferent E(B - V) values for the young and the old stellar population. 

^ ICN02 used a value of 68%. Here, conservatively, 50% is assumed, i.e., 
less of the ionizing radiation is reprocessed into the EBL. 
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Here, for simplicity, a single attenuation curve for all stellar ages 
specified by one E(B - V) value is used. Tiie resulting attenuation 
is in reasonable agreement with more detailed attenuation curves, 
e.g., the attenuation curves derived by |Pei| ( [T992j l for nearby galax- 
ies. Dust attenuation is strongest in the UV to O where the EBL 
limits considered here are not strongly constraining, and, therefore, 
the exact shape of the energy dependent dust attenuation is not im- 
portant for the analysis presented here. FigurefTltop panel illustrates 
the effect of dust attenuation on the resulting EBL SED (dotted vs. 
solid line). 

The attenuated photons are reemitted in the MIR and FIR. 
The dust reemission is modeled using the galaxy IR SED templates 
from |Chary & Elbazl ( |2001| l covering a range of IR luminosity from 
10** to 10" Lq. The IR SED depends on the total IR luminosity. To 
map the galaxy IR SED templates to the SSP spectra, the total dust 
attenuated luminosity L^^^^ from the SSP spectra is scaled by a fac- 
tor /;r (in Mq) and matched to a corresponding integrated galaxy IR 
luminosity Lf^'- from the IR SED templates: 



tGAL _ f ^ J SSP 



(4) 



The corresponding galaxy IR template is then scaled by a factor 
il fir and added to the SSP spectrum. For values of fir x L^J^J^ ex- 
ceeding the range provided by [Chary & Elbaz| J200T) the high- 
est/lowest galaxy IR template is scaled accordingly.!^ 

fr is a free parameter of the model. Given that the SSP spectra 
are normalized to the mass of the stellar popu lation, fr should b e 



2002 



of order of a galaxy mass (i.e. ~ 10 Mg; see Kneiske et al. 
The effect of different choices for /;, on the resulting EBL SED 
are illustrated in Figure[2] increasing values of /;, increase the ratio 
between the the FIR and MIR EBL. 

The parameters E{B - V) and fr can, therefore, be used to 
change the overall shape of the EBL SED, i.e., to shift photons 
from the UV/O/NIR to the MIR or FIR. Since the aim of the paper 
is to produce conservative upper limits, E{B- V) and /> are used to 
construct an EBL model which has minimal tension with the EBL 
limits used in the analysis, i.e., which produces an EBL SED fol- 
lowing closely the limit. This will ensure conservative upper limits. 

The resulting EBL SED for the optimum choice of parame- 
ters, E{B - V) = 0.15 and fr = 3 x 10', is shown in Figure |2] 
It follows very closely the upper limits from ME 12. For smaller 
or larger values of /> the resulting EBL SED exceeds the EBL 
upper limits in the MIR and FIR, respectively (Figure [2] Upper 
Panel). A higher/lower E(B - V) leads to overall more/less atten- 
uation and, thereby, decreases/increases the UV/O/NIR relative to 
the MIR and FIR emission (Figure |2] Lower Panel). The effect of 
different choices of E(B - V) on the limits will be discussed in 
Section|4] 



2.3 Star formation rate density 

The SFRD, as derived from a large compilation of different mea- 
surements by [Hopkins & Beacom,(2006, HB06 in the following), 
together with two recent estimates by Rodighiero et al.[ j2010| 
ROlO in the following) (IR) and |Karim et al.|p011| KAll in the 
following) (radio) is shown in Figure [3] The SFRD up to redshifts 
of z ~ 4 is characterized by a steep rise up to redshifts of z ~ 1 
followed by a flatter behavior up to redshifts of z = 3-4. While the 
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Figure 3. Star formation rate density (SFRD) vs redshift. Shown are the 
best fit broken power law function and the 1 and 3 <t contours from HB06 
and measurements from [Rodighiero et al.H20ro) and jKarim et al.|J2011^ . 
The dashed and dotted lines illustrate the behavior of a power law function 
in z + 1 for different choices for the power law index /?. For /? > the 
functions are normalized at z = 4 at the lower edge of the 3 o" contour from 
HB06, for /J < the functions are normalized at z = 0.9 to the measurement 
from , Karim et al.|j2011) . All values have been converted to a Chab IMF 
using a scaling factor of 1.65 (see Appendix [AJfor details). 



data agree on the overall shape, two different cases can be distin- 
guished, exemplarily represented by the ROlO and the KAl 1 data: 
in the former case, a higher peak normalization of the SFRD is ex- 
pected, followed by a flat or weakly declining SFRD for z > 1. 
In the latter case, the peak SFRD is lower (0(0.2 dex)) but then is 
followed by a still rising SFRD with a shallower slope up to red- 
shift z ~ 3 - 4 (a trend also found in H-ce data by [ Sobral et al.[ 
[2012[ >. At higher redshift the SFRD is less certain and different trac- 
ers produces different results: while observations of Lyman break 
galaxies in deep Hubble Space Telescope observations indicate a 
rapid decline for z > 4 (e.g. jBouwens et al.|2009 l SFRDs derived 
from gamma-ray burst observations indicate a flat behavior l [Yiiksel[ 
[et al.|2008[[Kistler et al.|2009| ). In this study it will be focused on 
the SFRD up to redshifts of z = 4, in particular on the peak of the 
SFRD around z~ 1. 

Overall, the SFRD p,(z) up to a redshift of z = 4 can be well 
described by a broken power law in z -i- 1 : 



pAz) = Po 



z+1 
Zn+ 1 



with 



a for z < Zo 
fi for z > Zo 



(5) 



(6) 



In the foffowing, values of/,- will be quoted without explicitly noting the 
unit. 



Zo being the redshift of the peak and po being the SFRD normaliza- 
tion at zo- 

This broken power law description of the SFRD is adopted for 
the EBL modeling. p,(0) is fixed to a value derived from the best fit 
broken power law from HB06: p,(0) = 0.02 and 0.012 for the Salp 
and the Chab IMF, respectively. Varying p,(0) in the 3 cr errors of 
HB06 (0(30 %)) does only result in small changes of the resulting 
EBL density (0(1 - 7 %), depending on the wavelength), mainly in 
the UV (and then, through dust attenuation and reemission in the 
MIR/FIR). Zo, Po, and/? are chosen as free parameters, which will 
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Figure 4. Limits on the SFRD parameters zo and po (ie- the peak red- 
shift and normalization) of the EBL modeling for /? = 0.3. The other 
model parameters are set to the fiducial values (Chab IMF, p,(0) = 
0.012Moyr-' Mpc^^ E{B - V) = 0.15, and /;> = 3 X lO'). Squares and 
crosses denote the values of zq and po for specific EBL models tested. z,q- 
po-pairs shown with a square symbol result in an EBL density in conflict 
with the upper limits (see legend in the lower right panel for details). 



be varied, while a is given by 



logPo - logp»(0) 
log(zo + 1) 



(7) 



The different values of/? adopted in this study are illustrated 
in FigurelS] Two different cases can be identified: a rising {/? < 0) 
and a falling (j8 > 0) second slope of the SFRD. For the former 
case, i.e., an increasing SFRD aX. z > Zo, P = -1 with a maximum 
redshift of Zmax = 3 is adopted, which is a good representation of 
the SFRD found by KAll for z > 0.9. In addition, /3 = -0.5 is 
used, which results in a reasonable representation of the SFRD for 
increased values of po at z ~ 1 above the value from KAll. For 
the latter case, i.e., a decreasing SFRD at z > zo, the choice for 
the different y8s is derived from the 3 cr contour of the broken power 
law fit to SFRD measurements from instantaneous tracers by HB06 
(Figure [3j: the lower edge of the 3 cr contour at z = 4 is used as 
normalization for the SFRD to derive different slopes which inter- 
sect the contours at z = 1 at the minimum (J5 = 0.3), the best fit 
(fi = 0.7), and the maximum (J3 = 1) allowed value. In addition, 
/3 = 1.5 is used which connects the minimum value at z = 4 with 
the maximum value at z = 2. As an extreme case not compatible 
with the SFRD values derived from instantaneous tracers j8 = 5 is 
considered, corresponding to a sharp drop of the SFRD after the 
peak. 



3 METHOD 

A grid in zo vs po is constructed to investigate different SFRDs in 
the EBL modeling, covering a range of z = 0.2 to z = 2.4 in redshift 
and log(p») = -1.2 to log(p,) = -0.5 in peak SFRD (see FigureHl. 
For each point (zo,Po) of the grid the EBL density is calculated. 
The resulting EBL SED is then compared with the EBL limits from 
ME 12 by calculating the ratio between the EBL SED P,,(0) and the 
limits La at different wavelengths: 



p.m 



(8) 



If the ratio t^ exceeds 1 at any wavelength, the EBL model 
parameters are considered to be in conflict with the limit. 

Two different cases can be discerned: (i) if the maximum ratio 
value exceeds one (i.e. max(f^) > 1) the resulting EBL is consid- 
ered in weak tension with the limits, (ii) Allowing for a systematic 
error of the EBL limit of 20% (see Appendix[B]for details), the EBL 
is considered to be in strong tension with the limits if the maximum 
ratio value at any wavelength exceed 1.2 {max(tA) > 1.2). 

Figure PI shows the results for a Chab IMF for/? = 0.3. The 
squares and crosses denote the value of zo andpo (i.e. the peak red- 
shift and normalization) for specific EBL model tested. The marker 
symbols and colors give the ratio of EBL SED vs limits (see the 
legend in the lower right panel of the Figure). Small squares repre- 
sent weak tension with the upper limit, while big squares represent 
strong tension. 

The constraints on zo and po are converted into constraints on 
the SFRD as follows: for each zo value tested the p^ value is cal- 
culated for which max(f,i) = 1/1.2, i.e., for which the weak/strong 
tension limiting value is reached. This is done via linear interpo- 
lation between the discrete values in t^ vs po. For each value pair 
(zo,Po ) the SFRD function p(z, Zo.Po ) is evaluated. The limit on 
the SFRD for a specific redshift z, is then taken to be the maximum 
of the different p{Zi,Zo,PQ ) functions, i.e., an upper envelope is 
calculated. 



4 RESULTS 

SFRD p, / f3 The upper limits on the SFRD for different values 
of /? are summarized in Figure B] First, it can be noted that for 
all choices of /3 compatible with the SFRD values from instanta- 
neous tracers the SFRD is constrained to < 0.1 Mg yr"' Mpc"' and 
< 0.2Moyr"' Mpc"' for z ~ 1 and z ~ 2, respectively. For z ~ 1 
this value is below the best fit value from HB06 for their broken 
power law fit. For p < 0.3 the 1 cr range from HB06 is in weak 
tension with the limits. The measurements from KAll are not in 
conflict with the upper limits derived here, even for a rising SFRD 
(/? < 0) beyond the peak redshift zo- In the following, /? = 0.3 will 
be adopted as the fiducial value, which is also the best fit value 
found by Ytiksel et al.| ( |2008| l in their broken power law fit of a re- 
cent SFRD data compilation including data from gamma-ray burst. 
Higher values of/? (> 1.5), which are in conflict with the values de- 
rived from instantaneous tracers (Figure [Sl, would lead to weaker 
limits. But even for the extreme case of/? = 5, corresponding to an 
almost instantaneous drop of the SFRD behind the peak, the limits 
are still in weak tension with the HB06 values. 

Metallicity Figure [6] displays the resulting upper limits on the 
SFRD for different metallicities of the SSR With decreasing metal- 
licity the constraints get significantly stronger, with the limit being 
~ 0.2 dex lower for Z = 10"* in comparison to the limit for the fidu- 
cial metallicity value of Z = 0.02 = Zq. This is a direct result of 
the increased UV, O, and NIR emission for the sub-solar metallic- 
ity models (FigurefTlbottom panel). Since the metallicity increased 
over the history of the universe, i.e., higher redshift should have 
lower metallicities, adopting solar metallicity as fiducial value is a 
conservative choice. For a metallicity higher then the solar value 
the limits are very similar to the solar metallicity ones, with only 
small difference at at low redshift z < 0.75, where the UV emission 
from young stars provides a larger contribution to the total EBL. 
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Figure 5. Limits on the SFRD for different values of /). Other EBL 
model parameters are fixed to the fiducial values (Chab IMF, p»(0) = 
0.012Moyr-' Mpc"^ £(B - V) = 0.15, and /, = 3 X lO'). The light 
red filled area coiTesponds to the exclusion region for the fiducial value 
/? = 0.3. The blue line and filled areas are the best fit and the 1 and 3 o" con- 
tours from HB06, respectively. Grey markers are the measurements from 
ROlO and KAl 1 (see Figure[5]for details). 




Figure 6. Limits on the SFRD for different metallicities of the SSP spectra. 
Other EBL model parameters are fixed to the fiducial values (Chab IMF, 
p,(0) = O.OllMoyr-' Mpc-3,/3 = 0.3, E{B-V) = 0.15, and/y = 3x lO'). 
The light red filled area corresponds to the exclusion region for the fiducial 
metallicity value Z = 0.02 = Z©. Other markers are as in Figurels] 
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Dust absorption The effect of different strengths of the dust atten- 
uation parameter E{B- V) on the resulting SFRD limits is displayed 
in Figure [7] Varying E{B - V) from the fiducial value by ±0.05 
does not strongly change the results: A larger value E{B - V) = 0.2 
leads to stronger constraints in the redshift range up to z ~ 1.5 and 
marginally relaxed limits at higher redshift. Lowering the value to 
£(5 - V) = 0.1 results in stronger constraints at redshift z ^ 0.5 
and marginally weaker limits below. This behavior can be under- 
stood from the correlation of the peak position of the EBL SED 
with the peak of the SFRD: a low-redshift peak of the SFRD leads 
to lower-wavelength UV/O peak of the EBL SED. While in the 
MIR/FIR the rising (low-wavelength) flank of the EBL SED peak 
is strongly constrained by the EBL limits the UV/O EBL is not 
strongly constrained. A model with a lower £(6- V) value produces 
less MIR/FIR EBL and more UV/O, which leads to a relaxation of 
the limits for lower-redshift-peaked SFRD. A similar argument can 
be made for stronger dust attenuation and higher SFRD peak red- 
shifts. For E{B - V) values ±0. 1 from the fiducial value the same 
trend are visible but the limit overall get stronger. 



Figure 7. Limits on the SFRD for different strengths of the dust attenuation 
parameter E{B - V). Other EBL model parameters are fixed to the fiducial 
values (Chab IMF p.(0) = 0.0I2Moyr-' Mpc"^ yS = 0.3, and /;> = 3 X 
lO'). The light red filled area corresponds to the exclusion region for the 
fiducial dust attenuation value E[B - V) = 0.15. Other markers are as in 
FigurelS] 



Initial mass function As discussed in Section lZTl a Chab IMF has 
been adopted as fiducial value for the EBL modeling, but indica- 
tions for variations of the IMF do exist. In FigurelSlthe SFRD limits 
for an EBL model with an Salp IMF are shown for different choices 
of /3. In comparison to the SFRD derived from instantaneous trac- 
ers the limits are more constraining then in the case of a Chab IMF 
effectively excluding the full 3 a range from HB06 at z ~ 1. 

This difference can be understood from the time evolution of 
the ratio of the SSP emissivities for a Salp and a Chab IMF: while 
for young ages of the SSP < 10' yrs the ratio for the luminosities 
between Chab and Salp IMF is /imf ~ 1.65 it decreases for larger 
ages (see Figure [AT] and discussion in Appendix [A). The EBL ef- 
fectively samples a large time range of order of the age of the uni- 
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Figure 8. Limits on tlie SFRD for a Salp IMF for different values of 
yS. Other EBL model parameters are fixed to fiducial values (p,(0) = 
0.02 Mo yr-' Mpc"\ E(B - V) = 0.15, and fir = 3 X lO'). The light red 
filled area corresponds to the exclusion region for the fiducial value /? = 0.3. 
Other markers are as in Figurep] 
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Figure 9. Ratio of the EBL SED for a Chab and Salp IMF Shovi^n are results 
for two different values of p.(0) for the Chab IMF: 0.012 Moyr"' Mpc"^ 
the fiducial value, and 0.02 Mq yr"' Mpc"^, the Salp IMF value. 



Figure 10. Number of EBL models excluded at a certain wavelength 
Nex('I) divided by the total number of excluded models Ng^^. In the up- 
per panel a Chab IMF in the lower panel a Salp IMF is used for the EBL 
modeling, respectively. The zo/po grid from Figureplis adopted and differ- 
ent values of E(B - V) are used. Other EBL model values are kept to the 



fiducial values {p.(0) = 0.012/0.02Moyr- 
/? = 0.3, and /,- = 3x lO'). 



Mpc-^ for a Salp/Chab IMF, 



verse. Since the SFRD displays a peak at a redshift of z ~ 1, the 
EBL at z = is already influenced by the emission from older stars, 
resulting in a decreased ratio for the EBL (see also |Horiuchi et af] 
|2009[ Section C). This effect can also be seen in Figure[9] where the 
ratio between EBL SEDs calculated with a Chab and a Salp IMF 
are shown. For typical model assumption for the SFRD the ratio 
in the EBL SED is between 1.3 and 1.4 in the wavelength range of 
1 - lOfiin, i.e., much lower then the 1.65 used as IMF scaling factor 
for the SFRD from instantaneous tracers. [Horiuchi et al.H2009l l. 

EBL wavelength dependency Since the EBL limits cover a range 
of over two decades in wavelength it is interesting to investigate at 
which wavelength the tension between the limits and the modeling 
actually occur. Several parameters entering the EBL modeling have 
no or only minor impact on the resulting EBL SED shape: chang- 
ing po changes the overall normalization of the resulting EBL (see 
Equation [T| and [2l. The effect of varying zo on the EBL SED shape 
is more complicated, mainly resulting in a shift of the EBL SED 
due to the increased redshift of the epoch with the maximum emis- 
sion. Given the limited range for zo considered here (~ 0.25 - 2.25, 
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i.e., a maximum shift in wavelength by a factor of three) the effect 
is also not overly strong but should be kept in mind. 

To investigate the EBL wavelength dependency the relative 
rejection at a specific EBL wavelength is calculated using a set 
of models marginalized over po and zo- The number of EBL mod- 
els excluded at a certain wavelength Nex(/1) is divided by the total 
number of excluded models Ng^^. Only models in strong tension 
are considered (though the results do not significantly change when 
also including models with weak tension). The results of this anal- 
ysis are shown in Figure [To] for a Chab (upper panel) and a Salp 
IMF (lower panel) and different values of E(B - V). 

It can be seen from the upper panel that for the fiducial choice 
of parameters (i.e. E(B-V) = 0.15) a larger fraction (~ 98 %) of the 
shapes are excluded in the NIR, MIR, and FIR simultaneously, i.e., 
the methods presented here take full advantage of the wide wave- 
length range of the EBL upper limits. This behavior is expected 
from the construction of the model: the NIR/MIR/FIR ratios have 
been tuned to closely follow the EBL upper limits. A lower value 
of E(B - V) leads to more exclusion in the NIR, a higher value to 
an increased exclusion in the FIR. The results using a Salp IMF 
are very similar, though a slightly higher value of E(B - V) some- 
where between 0.15 and 0.2 should lead to a more homogeneous 
exclusion over the wavelength range. The simultaneous tension in 
the three wavelength regimes (NIR/MIR/FIR) make the SFRD lim- 
its robust against, e.g., changes in the EBL upper limits caused by 
the exclusion of individual data or methods in the EBL upper limit 
calculations (see ME 12 for details). 



5 DISCUSSION AND CONCLUSIONS 

As discussed above, recent limits on the EBL density derived from 
HE/VHE spectra of distant sources constrain its density to a level 
close to the lower limits derived from integrated galaxy counts 
(ME 12). This implies that the bulk of the EBL density is produced 
by stellar emission (and stellar emission reprocessed by dust). The 
EBL density is, thus, directly connected to the star formation his- 
tory as determined by the SFRD. 

When converting these EBL limits into constraints on the 
SFRD via EBL modeling a peak SFRD lower then the one derived 
from instantaneous star formation tracers is implied. While the ten- 
sion between the SFRD from instantaneous tracers and the con- 
straints derived from the EBL limits for the fiducial model adopted 
in this study is not overly strong (~ I cr) it must be noted that the 
parameters of the EBL modeling have considerable uncertainties. 
For the study presented in this paper, they have been fixed to the 
most conservative values in the sense that they produce the weakest 
constraints (i.e. the lowest EBL densities). Consequently, the ten- 
sion will increase when the parameters are varied from their fiducial 
values, as it can be expected in the case of, e.g., the metallicity or 
possibly the IMF. 

One possibility to weaken this tension is to adopt an IMF 
which produces different luminosity outputs for the instantaneous 
and the late-time SFR tracers, which would imply a mid-heavy 
(1 - 8 Mo) IMF jFardal et al.|2007| >. Such an IMF could also ex- 
plain the difference between the SFRD predicted by instantaneous 
star formation tracers in comparison to the one derived from the 
evolution of the stellar mass, as traced by old stellar populations 
( [WiUdns et al.|[2008| . The main problem with such a mid-heavy 
IMF is that it is in contradiction with the finding of Ivan Dokkum 



( |2012^ , which both require additional contribution from low mass 
stars with masses < 0.5 M©. 

While, overall, the constraints on the EBL density are robust 
(see the discussion in ME 12), in the FIR wavelength regime they 
depend on a single measurement. When removing these data from 
the sample the EBL limits in the MIR to FIR are significantly weak- 
ened (see Figure 7d in ME12). Consequently, this would enable one 
to weaken the limits on the SFRD derived in this study by adopt- 
ing a stronger dust attenuation (a higher value of E(B - V), see 
Figure|2l. 

Other sources, e.g., active galactic nuclei or stars in the early 
universe, do also contribute to the overall EBL density. For exam- 
ple, [Dominguez et al.| ( |2011[ > find a contribution of 6 to 13 % to 
the bolometric EBL for AGN-type galaxies. Including such contri- 
butions in the modeling and, thereby, increasing the resulting EBL 
density for a given model, would further strengthen the limits. 

In this study, only upper limits on the SFRD have been de- 
rived. In principal, with a reasonably well constrained EBL den- 
sity through upper and lower limits, it is possible to derive a best 
fit SFRD. The main problem for such an ansatz is the insufficient 
knowledge about many of the parameters in the EBL modeling, 
e.g., the metallicity evolution, the dust attenuation etc. There is 
promising progress in reducing the uncertainties, e.g., from the ob- 
servations of gamma-ray burst afterglows which provide a direct 
view into the environment of stellar formation at higher redshifts 
(e.g. |Zafar et al.|201 l||Schady et al.|20I2[ ), but these measurements 
are often subject to strong observational biases and the extrapola- 
tion needed for a global modeling of the SF are still large. 

Upgrades of the current-generation VHE instruments 
(H.E.S.S. II, MAGIC-II) are being commissioned specifically 
targeting the sub- VHE (E < 100 GeV) regime, which is key for 
future VHE EBL studies. Currently under development, with 
deployment starting as early as 2014, is the Cherenkov Telescope 
Array (CTA, see |CTA Consortiumi 2010). CTA will deliver a 
ten-times improved sensitivity and extended energy coverage over 
current generation VHE instruments. EBL studies are one of the 
major science drivers for CTA and a precise measurement of the 
EBL density and its evolution from VHE observations will be in 
reach I jRaue & Mazin|2010^ . In the ultra-high gamma-ray energy 
regime (E > lOTeV) the non-imaging ground-based detector 
system SCORE is under development enabling a new view on the 
MIR to FIR EBL jTluczykont et al.|2011> . 



APPENDIX A: IMF CONVERSION FACTOR 

Assumptions about the IMF are necessary to convert measured 
quantities like luminosities into specific star formation rates and, 
thereby, the SFRD (see, e.g., HB06). For most quantities, like the 
IR or the radio luminosity, conversion factors have been calculated 
assuming a canonical Salp IMF ^Kennicutt| 1 998l l . As discussed in 
Section [2TT| the Chab IMF provides a reasonable description of the 
average IMF in the local universe. In the following, a scaling factor 



/iMF is derived to convert SFRDs calculated for a Salp IMF p, 
a Chab IMF pp"'': 



Salp 



to 



Chab _ Salp/ /■ 
P. - P* //IMF 



(Al) 



|& Conro}^(|2010|) and (maybe to a lesser degree) [Cappellari et al. 



For this, the specific luminosities ratios for different SSP spectra 
from BC03 for a Chab IMF and a Salp IMF are investigated. 

In the top panel of Figure [AT] the specific luminosities ratios 
from SSP spectra for a Chab and Salp IMF for different ages are 
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Figure Al. Upper Panel: SSP luminosity ratio for a Cliab and a Salp IMF 
for different ages of the stellar population. Lower Panel: SSP emissivity 
ratio for a Chab and a Salp IMF for different times assuming a constant 
SFR. 



shown. For young ages < 10^ yrs the ratio shows a slightly de- 
creasing slope from ~1.7 at ~ 0.05 yum to ~ 1.6 at ~ 20 fjm. For 
later ages the ratio decreases and shows a flat behavior in wave- 
length. In the lower panel the emissivity for different integration 
times are shown assuming a constant star formation rate. Here, the 
ratio shows a similar behavior as for the young age SSP spectra ra- 
tio, though with slightly lower values. In HB06 the emissivity value 
at O.ljxm is adopted as a good proxy for the overall scaling. Here 
this implies an IMF conversion factor of /imp = 1-65 (dotted line in 
the Figures). It can be seen from Figure [AT] that an IMF conversion 
factor of /imp = 1-65 is in good agreement with the young age SSP 
spectra and the constant star forming SSPs, though a wavelength 
dependent spread of about +5 % remains. 

[Horiuchi et al.|p009H derive a slightly higher conversion fac- 
tor of /imp - 1-8 for a Kroupa to a Salp IMF using the PEGASE.2 
code l |Fioc & Rocca-Volmerange|1997^ . 

For completeness, the ratio of the integrated ionizing luminos- 
ity and the luminosity absorbed by dust for the Chab and Salp IMF 
is shown in Figure \A2\ These quantities are directly related to the 
total infrared luminosity (i.e. in the EBL model the IR luminosity 
is the sum of the dust attenuated plus 50 % of the ionizing lumi- 



Figure A2. Ratio of the integrated ionizing luminosity and luminosity ab- 
sorbed by dust for SSP models with Chab and Salp IMF. 



nosity). Again, for ages up to 10' yrs, the ratio lies between 1.6 and 
1.7, i.e., /imp = 1.65 is a good choice. 



APPENDIX B: 
EBL LIMITS 



SYSTEMATIC UNCERTAINTIES IN THE 



The estimates for the systematic uncertainties in the EBL limits in 
ME12 are taken from |Mazin & Raue| ( |2007} , who used a similar 
technique as ME 12 to derive limits on the EBL density. The sys- 
tematic uncertainties are estimated to be ~ 32 % in the O to NIR 
and 33 - 55 % in the MIR to FIR. About ~ 30 % of this uncer- 
tainty results from the choice of the grid spacing, where a smaller 
grid spacing (i.e. allowing for sharper features in the EBL SED) 
leads to a relaxed limit. Additional errors arise from the neglected 
EBL evolution (~ 10%), and the uncertainties in the absolute en- 
ergy scale of ground-based VHE instruments taken to be ~ 15% 
(~ 10 - 45 % for the 0/NIR and the FIR/MIR, respectively). In 
the study of ME 12 the EBL evolution has been included, therefore 
this contribution can be neglected. In addition, recent studies show 
that the uncertainty in the absolute energy scale of ground-based 
VHE instruments is only at the level of ~ 5 % (Meyer et al. 2010|l. 
The width of the features in the EBL density produced by the EBL 
modeling in this paper are significantly larger then the ones used to 
derive the limits in ME 12 and, therefore, the systematic error aris- 
ing from a the choice of the grid spacing can safely be neglected in 
this study. It is, therefore, concluded that the conservative system- 
atic error on the EBL upper limits relevant for the study presented 
in this paper is at most ~ 20 %, resulting only from the remaining 
uncertainty in the absolute energy scale. 
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